Astronomy & Astrophysics manuscript no. AA' 18689 


©ESO 2012 


February 21, 2012 





Probing the envelopes of massive young stellar objects with 
diffraction limited mid-infrared imaging"^ 

H.E. Wheelwright^'^, W.J. de Wit^ R.D. Oudmaijer^, M.G. Hoare^, S.L. Lumsden^, T. Fujiyoshi^, and J.L. Close^ 

^ Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 

e-mail: hwheelwright@mpifr-bonn.mpg . de 
^ School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, UK 
^ European Southern Observatory, Alonso de Cordova 3107, Vitacura, Santiago, Chile 

^ Subaru Telescope, National Astronomical Observatory of Japan, 650 North Aohoku Place, Hilo, HI 96720, USA 



o 

<N 



<N 

' Received month dd, yyyy; accepted Month dd, yyyy 

^ . ABSTRACT 

QJ ' Context. Massive stars form whilst they are still embedded in dense envelopes. As a result, the roles of rotation, mass loss and accre- 

[T I ^ tion in massive star formation are not well understood. 

>A/Vt?s. This study evaluates the source of the Q-band, Ac - 19.5 jim, emission of massive young stellar objects (MYSOs). This al- 
lows us to determine the relative importance of rotation and outflow activity in shaping the circumstellar environments of MYSOs on 
1000 AU scales. 

Methods. We obtained diffraction limited mid-infrared images of a sample of 20 MYSOs using the VLT/VISIR and Subaru/COMICS 
instruments. For these 8 m class telescopes and the sample selected, the diffraction limit, ~0.6", corresponds to approximately 
1000 AU. We compare the images and the spectral energy distributions (SEDs) observed to a 2D, axis- symmetric dust radiative 
' transfer model that reproduces VLTI/MIDI observations of the MYSO W33A. We vary the inclination, mass infall rate, and outflow 

\ opening angle to simultaneously recreate the behaviour of the sample of MYSOs in the spatial and spectral domains. 

O I Results. The mid-IR emission of 70 percent of the MYSOs is spatially resolved. In the majority of cases, the spatial extent of their 

I ■ emission and their SEDs can be reproduced by the W33A model featuring an in-falling, rotating dusty envelope with outflow cavities. 

O ' There is independent evidence that most of the sources which are not fit by the model are associated with ultracompact H ii regions 

I and are thus more evolved. 

C/^ . Conclusions. We find that, in general, the diverse ~20 yt/m morphology of MYSOs can be attributed to warm dust in the walls of 

■ outflow cavities seen at different inclinations. This implies that the warm dust in the outflow cavity walls dominates the Q-band 

emission of MYSOs. In turn, this emphasises that outflows are an ubiquitous feature of massive star formation. 

Key words, stars: formation-stars: imaging-stars: early-type-stars: winds, outflows-infrared: stars-ISM: jets and outflows 
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, 1, Introduction suggestion that massive SF is intrinsically different to the low 
tH- . ^ ^ t , . t . , t mass case. However, recent radiation-hydrodynamical simula- 
^. Massive stars, M ^ 8 Mo, can play a dominant role in the evolu- ^j^^^^ ^j^^j demonstrate that the radiation is actually chan- 
cy tion of their host galaxies via feedback phenomena such as ion- ^^^^^^ ^^^^^^ i^^jl ^^^-^ c avities v^hile accretio n con- 
O izing radiation, strong stellar winds and supernova explosions, ^j^^^^^ ^^^^^^^ ^„ equatorial disc ('Yor ke & SonnhalT^ I2OO2L 
Cn^ . These processes exert a large influence on the outcome of sub- [Krumholz et al. 2009; Kuioer et al. 2010). Therefore, deviations 
7^. sequent star formation (SF). Massive SF can therefore be con- ^^^^ j^^^^^^j ^y^jnetry, which are sugge sted by the sp ectral 
> sidered as an integral part of the overall SF process and a key energy distribu tions of young massive stars dHarvev et all 19771: 
ingredient for a complete description of it. However, our com- louertler et al.l[T99l . remove the limit imposed by spherical ac- 
^ . prehension of many aspects of massive SF is still unsatisfactory. ^^^^-^^ Nonetheless, confirming that massive stars build up their 
d • ,f Tl^' " r" """"'T^ ^ final mass via disc accretion by direct observations is challeng- 
■ ■ ■ ''^tl^ ^f^'^""^ ^° generally much better understood process ■ ^^^^^^^ ^j^U^ ^^-^^^^ (MYSOs) are generally lo- 
of SF at low masses or not (see e.g. Zinnecker & Yorke 2007). ^^^^^ distances and their AU-scale circumstellar envi- 
The uncertainties regarding massive SF are in part the re- ronment remains unresolve d with single dish tele scopes at most 
suit of the short Kelvin-Helmholtz (KH) timescale associated wavelength regions (see e.g. lCesaroni et all2007l) . 

with massive stars: they contract rapidly to sustain a high lu- , , ,■ ■ r , ■ ^ ■ j- , j- 

u-i t-ii • • • Tf Long-baseline infrared interrerometry IS one of the few tech- 

minosity while still acquiring their mass. If the accretion pro- . , . , , •„• , . ^ 



, , . „ ^, . ^. 1 , T niques which offers the milli- arc second (mas) resolution re- 

ceeds spherically, the ensuing radiation pressure could limit the ^ , , ^ Tv/r^^o^^ i- / iT^~rr~ nr^TTTT] 

r; 1 . 11 . • ^ 1 /lA n P~c7 — quired to probe tor MYSO discs (see e.g. IFollert et al.l 120 IG: 

nna stp mr mass to annroYimatp v di) H .arson Starrri pn r-" t-i i 1 i ■■ — i i — ^ ' ' i " 



n 1 . • ^ Af\ A/f n P~c7 — quirea to probe tor MY^U aiscs (see e.g. iroiiert et ai.l izuiu 



ITOTII lKahjll974 Iwolflri-, ft Ca»rindlil ll98Tll, This kdl^Tlto I VW- f jilj I zyiU Pf w.| ;1 m i zyi Jnwilmjjn [ zyi I ) 

In par ticular, an exemplary result is delivered by IKraus et al.l 
* This paper is based on data obtained using the ESO VLT at the ^20ii using the Very Large Telescope Interferometer (VLTI) 
Paranal Observatory with programme 083.C-0795 and the Subaru tele- and the near-mfrared beam-combiner AMBER. These authors 
scope, which is operated by the National Astronomical Observatory of reconstruct a ^-band synthesis image of a 20Mq MYSO with 
Japan. an angular resolution of 2.4 mas. The resultant image shows hot 
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material in a disc-like geometry of approximately 20 AU in size. 
As a result, a consensus is beginning to emerge that the disc ac- 
cretion scenario can account for the formation of stars up to at 
least ~ 30 Mq. This implies that the MYSO circumstellar envi- 
ronment deviates significantly from spherical symmetry. In this 
scenario, rotation, accretion discs and outflow activity play a cru- 
cial role in providing the asymmetric environment required to 
facilitate accretion (see e.g. Beuther et al..i,2002) . 

The combination of a high stellar luminosity with accretion 
and outflow activity results in the circumstellar environment of 
MYSOs being shaped by several forces. Determining the rel- 
ative importance of these phenomena in sculpting the circum- 
stellar environments of MYSOs requires spatially resolved ob- 
servations. In particular, the mid-infrared (MIR) wavelength re- 
gion provides a wealth of information on the geometry of the 
accretion environment. Furthermore, this information can be ac- 
cessed by both interferometric techniques and diflraction limited 
single-dish imaging with 8m class telescopes. MIR interferomet- 
ric observations reveal that the A/^-band emission of MYSOs on 
scales of 50 t o 100 AU is dominated by warm dust located in 
the en velope (IdeWit et al.ll2007l: iLinz et al.l l2009l: IVehoff' et al.l 
I2OIOI) . It is expected that such envelopes will contain outflow 
cavities as there are many detections of outflow activity associ- 
ated with MYSOs (e.g. G35.20-0.74 in De Buizer 2006; IRAS 
20126-h41 04 in De Bu i zer 20 07; CepA HW2 in de Wit et al. 
2009). In Ide Wit et al.l (l2010l) . we argue that, specifically, the 
warm dust responsible for the A/^-band emission of W33A is lo- 
cated in the envelope close to the walls of the cavities evacuated 
by the outflow. 

The aim of this paper is twofold. We aim to assess the source 
of the MIR emission of a sample of MYSOs and thus test the 
premise that their MIR emission is dominated by warm dust 
located in outflow cavity walls. This will then allow us to de- 
termine whether outflows play an important role in shaping the 
environments of MYSOs (as opposed to, for example, rotation). 
We address these goals by comparing spatially resolved images 
of a sample (20 objects) to appropriate mo dels. Our approach is 
simila r to the work presented previously in lDW09: de Wit et~an 
(120091 hereafter DW09) with the addition of 2D modelling. We 
present MIR VLT images observed at 20yum for a sample of 
MYSOs. We then compare the spatially resolved images to two- 
dimensional, axis- symmetric dust radiative transfer models that 
feature rot ating, collapsing enve lopes and outflow cavities (the 
models of IWhitney et al.l [2003allbb . The paper is structured as 
follows. We describe the sample selection, the observations and 
data reduction procedures in Sect. [2l The MIR images are de- 
scribed in Sect.[3]and the radiative transfer model and its use are 
detailed in Sect. IH The results of the modelling are presented 
in Sect. [5] and discussed in Sect. [6l Finally, our conclusions are 
given in Sect. H 



2. Selection, observations and data reduction 

We aim to perform our analysis on a representative sample of 
MYSOs, in order to draw general conclusions. Selecting such a 
sam ple is not trivial . Initi al attempts employed IRAS data (see 
e^aUa & Stahledll99ll: iMolinari et al] 1 19961: ISridharan et al.l 
I2002h . As a result, such catalogues suff'ered from source con- 
fusion due to the large beam of IRAS . To rectify this problem 
and generate a Galaxy-wide, unbiased, MYSO sample, we have 
conduct ed a survey to detect and characterise MYSO s: the RMS 
survey (iLumsden et al. 1 120021: lUrauhart et ani2008h. Th e RMS 
is based on data of the MSX survev (lEgan et al.ll2003h . MSX 



data ofl'ers a significant improvement over IRAS data in terms 
of resolution (e.g. arcsecond rather than arcminute resolution), 
and thus enables the selection of a representative sample of MIR 
bright MYSOs. 

From the RMS database^ we selected objects with distance 
estimates within 3 kpc0 and with 21 yum MSX fluxes larger than 
30 Jy. The flux limit was imposed in order to obtain a decent 
signal-to-noise ratio in the wings of the resolved profiles. In to- 
tal, we observed 19 objects at 20 jim using the VLT Imager and 
Spectrometer for the Mid Infrared (VISIR, ^Laga ge et 



lager and 
al] l2004l) 



mounted at the Cassegrain focus of UT3 of the VLT (see Table[T]). 
Observations were conducted using the imaging mode of the in- 
strument and the Q3 filter which has a central wavelength of 
19.5 fim and a half-band-width of 0.4 jim. During the observa- 
tions, VISIR was equipped with a DRS 256x256 Si:As detector 
with an angular pixel size of 0.127''. This configuration enabled 
us to obtain oversampled, diffraction limited images, which with 
8 m class telescopes results in an angular resolution of approxi- 
mately 0.6''. 

In this paper, we also present the imaging observations and 
analysis of a key MYSO object W33A (G012.9090-00.2607). 
It was observed at the slightly longer wavelength of 24.5 jim 
with the CO MICS instrument (C ooled Mid Infrared Camera and 
Spectrometer lKataza et al.ll200Ql) mounted on the Cassegrain fo- 
cus of the Subaru telescope. We employed the imaging facility 
of COMICS and used a filter centred at 24.5 fim (see DW09 
for filter response functions). During the observations, COMICS 
was equipped with a 320x240 Si:As IBC detector which pro- 
vides oversampled diff'raction limited images with a pixel size 
of 0.13". A log of all the observations is presented in Table [T] 

For both the VISIR and COMICS observations, standard 
stars were observed to provide a reference point- spread-function 
(PSF). The standard stars selected are MIR bright, -100 Jy 
at 25 yum , isolated sources which are expected to be unre- 
solved. Data reduction, consisting of shifting and adding nod- 
ded and chopped images, was conducted using routines writ- 
ten in iDL. The resultant images were not astrometrically cor- 
rected. Therefore, the images are presented in terms of distance 
from the peak of the intensity distribution of the target source. 
The azimuthally averaged intensity profiles of the standards are 
displayed in Fig. [T] The radial profiles of the VISIR standard 
stars are generally consistent until a level of approximately 0.001 
times the peak ffux and appear to behave as point sources. The 
standard star of the COMICS data (HD 124897) is of lower 
signal-to-noise ratio but is similar out to 1 percent of the peak 
ffux to the standards presented in DW09. 



3. Observational results 

We assessed the MIR images of our target sources by comparing 
their azimuthally averaged intensity proffles to the instrumental 
PSF determined from the standard stars. We ffnd that the ma- 
jority of the MYSOs observed, 14 of 20, are spatially resolved. 
The images of the target MYSOs are presented in Figs. [2] and [3l 
which contain the unresolved and resolved MYSOs respectively. 

The resolved MYSOs generally appear as a single source 
within the VISIR ffeld of view, 30" x 30". Most display an 
approximately circular symmetric morphology at the resolution 



^ http://www.ast.leeds.ac.uk/RMS/ 

^ Some of distances estimates have since been revised to larger val- 
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Table 1. The log of the observations. 
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08 48 48 64 


-43 32 29.1 


09/05/1 1 

yJ y 1 yJ^l X X 


760 


1.1 


900 
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4.1 


93 


G263 7759-00 4281 


08 46 34 84 

VJ \J ^ VJ ^) ^ . \J ^ 


-43 54 29.9 


09/05/1 1+05/12 


1500 


1.1 


150 


1.5 


4.0 


65 


G265 1438+01 4548 


08 59 27 40 


-43 4i 03 7 


09/06/04+05/05 


1833 


1.4 


150 


2.5 


8.6 


42 


G268 3957-00 4842 


09 03 25 08 

vj y \j ^ .^■j.vjvj 


-47 28 27.6 


09/06/05+06/22 

\J y 1 VJVJ/ VJ,J 1^ VJVJ/ 


2x 1305 


1.6/1.7 


250 


1.6 


3.7 


50 


G269 1 586-01 1 383 


09 03 31 76 


-48 28 45 5 


09/05/1 2 


700 


1.2 


25 


2.7 


7.3 


98 


G310 0135+00 3892 


13 51 37 85 


-61 39 07 5 

\J± •JV VJ/.^ 


09/05/26 

\j y 1 \j ^ 1 VJ 


120 


1.4 


250 


3.3 


57 


259 
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14 26 26.28 


-60 38 31 5 

VJVJ ^7 *J ^ 1 . ^ 


09/05/26 


2000 


1.3 


600 


3.7 


13 


58 


G318 0489+00 0854 


14 53 42 99 

1 ^ ^^,yy 


-59 08 56 5 

VJVJ ^ VJ . ^ 


09/06/01+06/07 

\J y 1 \J\Jf VJ X VJ / 


2x 1996 


1.2/1.2 


50 


7 


? 


41 


G318 9480-00 1969 

VJ ^) V \j , y 1 u \j VJ VJ . i y \j y 


15 00 55 31 

1 ^ VJVJ w'w'.^JX 


-58 58 52 6 

^ VJ ■JO w'.^.xJ 


09/05/10 

y 1 1 X VJ 


2000 


1.2 


650 


2.6 


11 


55 


G326 4755+00 6947 

VJ,J^\J.T^ / ^ ^ \^ \J\J ,\J y \ I 


1 5 43 1 8 97 


-54 07 35 6 

^ 1 VJ / ^.J.W 


09/06/07 

VJ V/ \J\Jj VJ / 
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650 
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8.8 


42 


G332 l^M^^l 2799 

VJ ^ ^) ^ , ^ y 1 1 1 vj^ , ^ 1 y y 


16 05 41 84 


-49 11 29.5 


09/06/07 

y 1 yjyj/ VJ / 
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1.2 
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1.9 
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09/06/07 

VJ V/ \J\Jf VJ / 
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3.5 
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-50 39 12.6 


09/06/07 
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1.4 


50 


3.5 


2.9 


48 


G339.6221-00.1209 


16 46 05.99 


-45 36 43.9 


09/06/07 


2225 


1.4 


150 


13.1 


520 


39 


G341. 128 1-00.3466 


16 52 33.19 


-44 36 10.8 


09/06/08 


2423 


1.1 


400 


3.5 


6.2 


37 


G343.5024-00.0145 


16 59 20.90 


-42 32 38.4 


09/06/07 


400 


1.4 


25 


3.0 


18 


131 


G343.5213-00.5171 


17 01 34.04 


-42 50 19.7 


09/05/10 


1497 


1.1 


75 


3.0 


13 


47 


G345.0061+01.7944 


16 56 46.37 


-40 14 26.7 


09/05/10 


280 


1.0 


100 


1.8 


8.7 


154 


G349.7215+00.1203 
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09/06/07 
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1.6 


200 


24.4 


310 


31 


COMICS 
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09/08/30 
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14 06 40.95 


-36 22 11.8 


09/06/23 


2500 


1.0 


350 
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16 48 39.90 


-69 01 39.8 


09/05/10+05/26 


3600/3600 


1.4/1.4 


1800/600 






HD 163376 


17 57 47.80 


-41 42 58.7 


09/05/26 


3600 


1.1 


150 








HD 81797 


09 27 35.24 


-08 39 31.0 


09/05/12 


2500 


1.2 


850 








COMICS 




















HD 124897 


14 15 39.67 


+ 19 10 56.7 


09/08/30 


856 


1.9 


50 









Notes. The signal to noise ratio is defined as the ratio of the central peak to the root-mean- squared (rms) background noise. The distances and 
luminosities are taken from the RMS database. 
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Fig. 1. The azimuthally averaged intensity profiles of the PSF 
standard stars. HD 124897 is the PSF standard taken with 
Subaru/COMICS. The solid line with square points is the av- 
erage VISIR PSF. 



of the observations (~0.6'0, although some show clear signs of 
structure (e.g G263. 7759-00.4281, Fig. [3 panel A). Several ob- 
jects are associated with additional, distinct sources of localised 
emission, which often have a cometary type morphology (e.g. 
G269.1586-01.1385 & G349.7215+00.120 3, Fi^M panels D 
& M). This is suggestive of a Hii region (iHoare et al ][2007|). 



Indeed, both the objects mentioned above are detected at radio 
wavelengths at coordinates c onsistent with the loc ations of the 
cometary 20 jim sources (see lUrquhart et al.ll2007l) . It is of in- 
terest to note that 24.5 yum images of regions containing both 
MYSOs and ultracompact Hii regions show a clear morpholog- 
ical difiference between the two types of source, i.e. compact vs 
extended (DW09). 



Approximately one third of the sample are unresolved. Using 
the Kolmogorov-Smirnov test, we investigated whether the re- 
solved and unresolved samples were drawn from intrinsically 
diflTerent distributions of MYSOs in terms of distance and lumi- 
nosity. There is no significant diflTerence between the distance 
and luminosity distributions of the two samples, as might be ex- 
pected given that the objects are drawn from a single catalogue. 
However, the angular size of a centrally heated source of flux is 
dependent upon both the luminosity of central object and the dis- 
tance to it (angular size oc see e.g. lVinkovic & Jurkidll2007l) . 
We find that the resolved objects generally have a higher value of 

than the unresolved objects. More specifically, the hypothe- 
sis that the unresolved sources have the same distribution as 
the resolved sources can be discarded at a level of 98 percent sig- 
nificance. Thus, we conclude that resolved sources appear larger 
than the unresolved sources as the resolved sources are gener- 
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ally more luminous and less distant than their unresolved coun- 
terparts. 



4. Radiative transfer modelling 

To investigate the nature of the circumstellar structure of the re- 
solved sources, we compare the observed images to 2D axis- 
symmetric radiative transfer models appropriate for massive 
stellar objects surrounded by in-falling, rotating envelopes. By 
simultaneously comparing the resolved MIR images and SEDs 
to the radiative transfer calculations, we aim to understand the 
source of the 20yum MIR emission of MYSOs. The models and 
the modelling procedure are described in detail in the following 
sections. 



4.1. The radiation transfer code 

We employ the 2D, axis- symmetr ic dust radiation tran sfer code 
of Whitney et al. (for details see IWhitney et al.ir2Q03 a.b). The 
code calculates radiation transfer through a dusty structure 
which consists of a proto- stellar envelope surrounding a central 
star. A low density polar cavity can be inserted into the overall 
envelope structure. The code also allows for a dusty circumstel- 
lar disc. The density distribution of the envelope i s described by 
the treatment of a c ollap sing, rotating envelope of'U lrichI (Il976k 
and iTerebey et al.l (Il984i) . The key parameters determining the 
envelope density are the mass infall rate and the centrifugal ra- 
dius. The dust sublimation radius is calculated self-consistently 
from the stellar luminosity and the dust sublimation tempera- 
ture. The shape of the polar cavities can be described in two 
ways: using a polynomial function or following streamlines. The 
streamline is conical on large scales, while the apparent opening 
angle of the polynomial function can be specified (the opening 
angle at the stellar surface is 180°). If inc luded, the dust disc 
is flar ed and follows the a prescription (see lShakura & SunyaevI 
1 19731) to account for flux due to accretion. The reader is referred 



to 



Whitnev et all (l2Q03allbb for more details on the code and the 



possible geometries. We have used this code extensively and the 
reader is referred to de Wit et al. (2010, 2011) for more details 
on the application of the code to model MYSOs and their envi- 
ronments. 



4.2. Modelling methodology 

Our previous use of the Whitney et al. code consisted of con- 
structing a series of dedicated models in order to reproduce 
the SEDs, interferometric and auxiliary spatial information of 
MYSOs. We will exploit these customised models appropriate 
for MYSO environments in our approach here, partially mo- 
tivate d by the lack of suc h models in the well-known SED 
grid (iRobitaille et al.ll2007l) . In particular, we will use a model 
geometry that pro vides a successful fit to the source W33A 
(Ide Wit et~ani2010h . Its defining feature is that the MIR emission 
in the A^-band originates in warm dust located in the walls of out- 
flow cavities. The model proto- stellar envelope extends inwards 
down to the dust sublimation radius and it has low density out- 
flow cav ities as observed i n several MIR observations of MYSOs 
(see e.g. lDe Buizei]|2007b . We note that the model does not ex- 
plicitely include a dust disc (according to the prescription of the 
Whitney et al. code). This is motivated by the aforementioned 
MIR interferometric observations. Adding a disc generally re- 
sults in high visibilities in the MIR whereas MYSOs typically 
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Fig. 4. Panel a: Logarithmically scaled 24.5 yum image of the 
standard W33A model. Panel b: The same image as panel a) but 
convolved with the COMICS instrumental PSF. Panel c: The 
COMICS image of W33A. Panel d: The azimuthally averaged 
radial intensity profile at 24.5 yum of W33A (points with error 
bars) alongside the model radial profile and that of a PSF stan- 
dard (the lower line). The outflow axis of the model is aligned 
with the SE/NW orientation of W33As outflow. 



exhibit low visi bilities, indicative o f a more extended dust distri- 
bution (see e.g. lde WitetaDl2010l) . 

We illustrate the model 24.5 jim image for the W33A model 
in panel a) of Fig.|4j Also shown are the PSF convolved image 
(panel b), the COMICS image of W33A and the corresponding 
radial intensity profile, panels c) and d) respectively. The shape 
of the outflow cavities can clearly be delineated in the model 
image in panel a). This demonstrates that the dust located in 
the outflow cavity walls is warmed-up and emits a significant 
fraction of the MIR flux. W33A is know to d rive an outflow 
in the SE/NW direction (see iDavies et"aL and references 

therein). The orientation of the outflow is traced by nebulosity 
in the NIR. While the 24.5 yum image of this object is resolved 
(panel c), it appears relatively symmetric. However, the image 
is slightly extended towards the SE, in the direction of the blue 
lobe of its outflow. In panel d) we compare the model's radial in- 
tensity profile with the observed one and find they are essentially 
identical. We underline that we did not perform a model fit to the 
24.5 yum image, we simply used the final W33A model presented 
in de Wit et al. (2010) to create the corresponding image. 

Motivated by the good match of the W33A model and im- 
age and the knowledge that the SEDs of MYSOs are known to 
be similar, we proceeded to fit the observations of the remain- 
ing MYSOs in the following way. The SEDs of the MYSOs ob- 
served were constructed from the continuum data assembled on 
the RMS database (which mak es use of the work of various au- 
thors) and the methodology of iMottram et al.l (l201ll) . The SED 
coverage diff'ers from one source to the other. Generally one can 
distinguish between sources with continuum flux measurements 
up to 100 fim and those with sub-mm observations. We note that 
we do not fit explicitly the 9.7 yum silicate feature, contrary to the 
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modelling efforts in DW09. This is simply because this informa- 
tion is not available for the majority of our sources. 

To reproduce both the SEDs and 20yum intensity distribu- 
tions of the resolved sample, we used the W33A model geom- 
etry summarised above. We only altered parameters that can be 
expected to vary between objects: the inclination, the envelope 
infall rate (which scales the total dust mass in the envelope) and 
the opening angle of the outflow cavity (which may vary with 
time). We emphasise that we use one basic model to reproduce 
the observations. There may be differences between sources, for 
example the centrifugal radius and the properties of circumstel- 
lar discs - if any are present - may vary from source to source. 
The model will not account for such differences. However, it pro- 
vides a simple test of the hypothesis that the bulk of the Q-band 
emission of MYSOs can be attributed to warm dust in outflow 
cavity walls. Here we outline the methodology followed in re- 
producing the observations. 

For a given MYSO, the stellar luminosity was initially set 
to the value in Table [T] and the system inclination (the angle be- 
tween the polar axis and the line of sight) was estimated from 
its MIR image. The stellar luminosity and infall rate were then 
varied until the model SED was consistent with that observed. 
This set the luminosity and provided first estimates of the other 
parameters. Once the SED was satisfactorily reproduced, the in- 
clination and outflow opening angle were varied in an attempt to 
recreate the observed spatial intensity profile. This had an affect 
on the reproduction of the SED. Therefore, the infall rate was 
allowed to vary to enable a fit to both the SED and the radial 
profile. The model fitting was done by hand rather than using a 

minimisation technique as the computational time required to 
make a grid of model images for each object was prohibitively 
large. 

We demonstrate the effect of varying the free parameters in 
Fig. [5] In general, the more inclined models appear more ex- 
tended. This might be expected if much of the MIR emission 
originates in outflow cavity walls. At large inclinations, the pro- 
jected area of the outflow cavity walls is larger and thus the resul- 
tant flux distribution is more extended. The cavity opening angle 
has a similar effect on the intensity distribution. This is because 
increasing the opening angle increases the distance at which the 
majority of the cavity wall is exposed to direct irradiation from 
the central object. This also increases the projected area of the 
outflow cavities and the extension seen in the azimuthally aver- 
aged intensity distribution. Since both parameters affect the ob- 
served extension, they are slightly degenerate. However, the two 
parameters have different effects on the gradient of the radial 
profile (see Fig.O. 

The final SEDs are presented in Fig.[6]and the associated az- 
imuthally averaged intensity profiles are displayed in Fig. [71 The 
parameters of the selected models are presented in Table[2l Since 
the sample exhibits a varied morphology (see Fig.O, we discuss 
the modelling results for each object individually in Appendix 
IaI We discuss the general results in the following section. 



5. Modelling results 

In general, the observed SEDs could be reproduced relatively 
well. There are several cases where the observed NIR fluxes 
could not be reproduced exactly. However, we do not weight the 
NIR fluxes heavily as these can be strongly dependent on dust 
opacities and the circumstellar geometry on small scales. In the 
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Fig. 5. An exploration of the parameter space immediately sur- 
rounding the model of G310.0135-hOO.3892. The plots show the 
effect of varying the free parameters (the inclination, opening 
angle and infall rate). For comparison, a PSF standard is also 
shown in the intensity distribution figures. 



majority of cases, barring objects associated with Hii regions, 
the spatial distribution of the 20yum emission of the MYSOs 
observed can also be reproduced by our 2D axis- symmetric ra- 
diative transfer (RT) model. The individual model fitting re- 
sults for the SEDs and the 20 yum intensity profiles are shown 
in Figs.[6]and|7l Clearly, because of the unequal SED coverage 
from source to source, the model fidelity differs. We do not at- 
tempt to provide a unique model for each single source but eval- 
uate whether the W33A model is capable of providing satisfac- 
tory fits, by changing a few of the model's most basic parame- 
ters. From this starting point, we find that in general, if the image 
exhibits a single source with an approximately symmetric mor- 
phology, both the object's SED and the 20yum intensity profile 
are well reproduced. By extension of the case for W33A, the 
majority of the 20yum emission in these instances is very likely 
emission from warm dust in cavity walls. 

When an object's morphology appears either very complex 
or cometary, the model generally fails to reproduce the intensity 
profile. For example, this is the case for G343.5024-00.0145, 
which is a known ultracompact Hii region. For this object, a de- 
cent fit to the SED (panels I in Figs. [6]and[7]) predicts a "MYSO" 
structure that is much more compact than the one observed, even 
with extreme values for the inclination and opening angle. A 
cometary appearance alone will often indicate an ionized nature 
for the region, and this is generally confirmed via radio contin- 
uum observations. That the model cannot reproduce the obser- 
vations of MYSOs associated with Hii regions is perhaps to be 
expected as Hii regions represent a source of flux not considered 
by the model. Furthermore, the appearance of a Hii region likely 
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Table 2. The key parameters of the individual models. 
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Notes. We note that the output luminosity is affected by inclination 
angle, and is thus not necessarily identical to the source luminosity. 
Parameters not listed were not varied. In the case of several objects, 
particularly those associated with Hii regions, e.g. G343. 5024-00.0145, 
the model only reproduces the objects' SED. ^op.ang. denotes half the full 
opening angle. 



signifies the end of the embedded accretion phase of a MYSO, 
and it is this phase the model is designed to represent. 

We give an overview of the values for the varied parame- 
ters in TableO The opening angles of the best-fitting models are 
generally narrow (~ 10°). This implies that the outflows of the 
MYSOs observed are relatively well collimated. The envelope 
infall rates of the models, a parametrisation of the density distri- 
bution, are generally of the order lO'^^Moyr"^. The correspond- 
ing values of visual extinction. Ay, are typically bet ween 10 and 
30, in agreement with the measured AyS of MYSOs (iPorter et alJ 
Il998l) . although several of the more inclined models have Ay s in 
excess of 100. The total masses of the envelopes are generally 
of the order 10,000 M©. We note that these are dependent on the 
envelope outer radii, which are set at 5 x 10^ AU. In the case of 
the W33A model, reducing the outer radius by a factor of 2 de- 
creases the total mass by a factor of approximately 5. This is a re- 
sult of adopting a low but constant density, p = 8 x 10"^^g cm"^, 
within the outflow cavity. The envelope masses calculated are 
comparable to the c louds that the objects are located in (see 
lUrquhart et al.ll20lT[ and references therein). Even allowing for 
an uncertainty of a factor of a few, this indicates that these ob- 
jects are still heavily embedded in their natal material and their 
environment still contains a large reservoir of material available 
to form further stars. 



6. Discussion 

6.1. On the mid-IR emission of MYSOs 

Over the past decade, observational evidence indicating that 
the formation of massive stars, up to a mass of at least ~ 
30 Mq, is accompanied by phenomena characteristic of low 
mass SF has accumula ted. Such phenoni e na include: c i rcum- 
stellar discs (see e.g. Kraus et alj I201QI: iDavies et"al 
iMasque et al.l201 ll:[God di et al. 2011 ^, molecul ar outflows (see 
e.g. lBeut her et a l.ll2002|: [Cvganowski et al. 2009) a nd jet a ctivity 
(see e.g. lRodriguez & Bastianiil994i:iGuzman et ani201Ql) . 



Despite the many common characteristics of low and high 
mass SF, some observational phenomena seem to be exclu- 
sively associated with young, embedded high-mass stars. One 
example is the Class II 6.7-GHz methanol maser, whose ac- 
tivit y is uniquely assoc iated with massive star forming regions 
(e.g. lWalsh et al.|[T997l) . These sites evidently produce sufficient 
IR radiation in order to pump this transition, unlike sites of 
low mass SF. It was initially thought that this maser activity 
originated in circumstellar discs. As a result, a correlation be- 
tween the extension of several MYSOs in the MIR and their 
maser emission led to the sug gestion that the MIR e mission of 
MYSOs traces discs (see e.g. IPe Buizer et al.ll200Qh . However, 
observations with 8m class telescopes were able to spatially 
resolve the MIR emission of a number of masing sources to 
reveal that their MIR emission is aligned with, rather than 
perpendicu lar to, their large CO molecular outflows (see e.g. 
lDeBuizedl2006[l2007h . These findings associate the MIR emis- 
sion of massive YSOs with outflow cavities, rather than circum- 
stellar discs. Nevertheless, the kinematic diversity of methanol 
masers in a large sample of massive outflow sources does not 
support a uni quely identifiable source within the close MYSO 
environment (ICvganowski et al.ll2009t iMoscadelH et al.l l201 ih . 
Therefore, the source of the MIR emission of MYSOs has yet to 
be conclusively established. 

Recently, several massive young stellar objects with compact 
MIR emission were investigated using long-baseline interferom- 
etry. In several cases, the spatially resolved A/^-band emission on 
scales of 100 AU was found to be consistent with thermal enve- 
lope emission (de Wit et al. 2007; Linz et al. 2009) and in par- 
ticular with emission arising from the cavity walls (de Wit et al. 
2010). However, in some other cases, the mas- scale MIR emis- 
sion of MYSOs is suggest ed to be (at least partially) located in a 
circumstellar disc (see e.g. lFollert et al.l201Ql:lde Wit et al.ll201 ll: 
iGrellmann et al.l201 ll) . In other words, there are hints that the N- 
band emission, resolved with interferometry, is not completely 
dominated by envelope emission. The analysis of the Q-band 
emission presented in this paper favours the cavity wall emis- 
sion interpretation for a sample of MYSOs. Upon first inspec- 
tion, this appears at odds with the scenario of a disc dominating 
the N-band emission. Here, we consider the possible contribu- 
tion to the 2-band emission by a circumstellar disc in order to 
assess whether our finding agrees with theoretical expectations. 

The recent study by lZhang & TanI (1201 ll) discusses in detail 
the appearance of embedded young massive stars in the IR wave- 
length region. The authors use a dust radiative transfer code, 
as we do. However, their geometry differs from the one em- 
ployed in this paper in that their code includes an accretion disc 
that extends from the stellar surface to the centrifugal radius. 
Therefore, their results allow us to assess the possible disc con- 
tributio n to the Q-band emi ssion of MYSOs. The RT calcula- 
tions of IZhang & TanI (1201 ll) reveal that outffow cavities have a 
large inffuence on the shape of the SED. Furthermore, in their 
calculation for a 60^ inclination (see e.g. their Fig. 9), the out- 
ffow cavities are clearly signiffcant features in the images, up to 
a wavelength of 70 yum. In the NIR, scattered light from the cav- 
ities makes a signiffcant contribution to the near-IR emission. In 
the A/^-band, thermal emission from the disc and cental object are 
also evident. In the 2-band however, the emission from the base 
of the cavity clearly dominates the ffux. 

The relative contributions of the disc and envelope to the 
modelled SED depend critically on exactly those parameters 
with which we fft our observations, i.e. opening angle, mass 
infall rate and inclination. We note that the opening angle ap- 
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plied bv lZhang&T^ (l201lh is rather wide (2a ^ 90^), which 
would imply an object in a rather advanced phase of forma- 
tion, following the expec ted widening of the opening angle as 
function of time (see e.g. iKuiper et al.|[2blQl) . Smaller opening 
angles increase the extinction towards the central regions and 
thus favour the dominance of cavity wall emission in the total 
MIR emission. A similar reasoning applies for the mass infall 
rate. It is therefore not surprising that the observed MYSOs, 
which are believed to be in a less evo lved phase than the p hase 
represented by the fiducial model of IZhang & TanI (l201ll) . are 
dominated by the cavity wall emission. This argues against a 
large eff'ect, if any, by a circumstellar disc on our VISIR im- 
ages. This conclusion can also be inferred from the spectral en- 
ergy distribution of an accretion disc, since it peaks at wave- 
lengths much shorter than 20yum. Furthermore, the envelope 
contains much more cool material than does any reasonable ac- 
cretion disc, whose s izes for MYSOs are estimated to be less 
than 500 AU (see e.g. IPatel et al.ll2005l: lHoardl2006l: iReid et al 



20071: 



Jimenez-Serra et al.ll2007l:lde Wit et al.ll201lUGoddi et al 



2011bh . We conclude that it is reasonable to expect that dust 
emission from the envelope (the cavity walls) dominates the con- 
tinuum SED at wavelengths longward of the N-band. This is in 
agreement with our finding that the observed extension of the 
2-band flux can be modelled as emission from outflow cavity 
walls. 



6.2. The structure of MYSO envelopes: rotation versus 
outflows 

We have shown that the morphology of the 20 yum emission of 
MYSOs is consistent with that predicted by models incorpo- 
rating outflow cavities. We followed a similar methodology to 
DW09 and have improved on the modelling part of their analy- 
sis. DW09 compared their images to 1-D radiative transfer mod- 
els. This was motivated by the rather spherical appearance of 
most sources, which suggests that their MIR emission originates 
from a relatively symmetrical dusty envelope. They found that, 
in general, their data were best fit by an envelope with a density 
gradient of p oc r"^. This density law is shallower than that ob- 
tained via a similar analysis performed in the sub-mm (see e.g. 
iMueller et al.ll2002b . It was therefore suggested that this flatten- 
ing of the density law towards smaller spatial scales could be 
evidence for the onset of rotational support of the envelope at 
approximately 1000 AU. This suggestion is probably invalidated 
by the results presented in this paper based on more realistic RT 
models. 

Thermal emission from an envelope with evacuated outflow 
cavities fits the SED and the A^-band (de Wit et al. 2010) and Q- 
band (this paper) morphology of the luminous MYSO W33A. 
Additionally, in this particular case, the envelope model also fits 
the observed morphologies of the near-IR scattering nebula and 
the 350yum emission. This provides strong evidence in favour of 
this model and our explanation of the 20 yum emission (warm 
dust in the cavity walls). Moreover, the W33A image (Fig.|4l) is 
slightly extended towards the SE which is the direction of the 
blue lobe of the large scale outflow, reinforcing this interpre- 
tation. Importantly, we note that W33A is not the only case in 
our sample for which the extension at 20 jim is in the direc- 
tion of the outflow position angle (see for example the sources 
G263.7759-00.4281 and G332.2941-h02.2799). Therefore, our 
association of the observed MIR flux with outflow cavity walls 
can be substantiated. 



The brightness of the cavity walls, as exemplified in panel 
a) of Fig.m is a consequence of the density contrast between 
the in-falling envelope and the outflow cavities. The radiation of 
the central source cannot penetrate far in the dense envelope, but 
can heat dust in the walls of the outflow cavity to large radii. 
Once the model images are convolved with an instrumental PSF, 
it is no longer obvious that the flux traces outflows, despite the 
source being resolved. The flux from the outflow cavity walls 
can appear more extended than symmetric envelope emission. A 
centrally concentrated density distribution (steep density law) in 
a spherical envelope produces a small thermal emission region 
in the MIR. The shallower the density law the larger the (nor- 
malised) emission region becomes. Therefore, the influence of 
the outflow cavities can mimic the efl'ect of a shallow density 
gradient. We verified that the sizes of the emission regions dif- 
fer strongly between a -2 and -1 density law, but the diff'erence 
between a flat density law and a - 1 density law are relatively mi- 
nor. Still the SEDs difl'er significantly. The steeper density law 
produces higher MIR fluxes. As an exercise, we fitted the stan- 
dard 2D W33A model with a spherical model and find that the 
SED and intensity profiles are matched best by a density law be- 
tween -0 and -1. This explains the results presented in DW09. 
Given the more realistic nature of the 2D models we have em- 
ployed, we suggest that the preference for shallow density laws 
found by DW09 is the result of much of the MIR emission of 
MYSO emanating from warm dust in the cavity walls. An ex- 
planation in terms of rotational support on scales of -1000 AU 
cannot be substantiated. 



6.3. SED fitting 

Fitting spectral energy distributions has long been used as a 
tool to provide important constra ints on the morpholo gy of cir- 
cumstellar environments (see e.g. lGuertler et al.|[l99ll) . Still, the 
inclusion of spatially resolved observations at diff'erent wave- 
length regimes is indispensabl e for a co rrect interpretation of 
the observed emission (see e.g. Ivan der T ak et al. 2000). This is 
because SED fitting can result in degenerate solutions. To as- 
sess the added value provided by the VISIR im ages presented 
in this paper, we use the results presented in iMottram et al.l 
(I2OIII). These authors used the grid of models provided by 



iRobitaiUe et al.l (l2007t) to fit the SEDs of a large percentage of 
MYSOs from the RMS survey in order to determine t heir bolo- 
metric luminosity. The grid of iRobitaille et al.l (l2007l) was cre- 
ated using the same RT code of Whitney et al. which we use 
in this paper. Therefore, we us ed the code to create images for 
the best fitting models found in lMottram et al.l (1201 ll) (in the fil- 
ter used) and assessed whether the model images are consistent 
with our observations. 

To illustrate the results of this exercise, we focus on the ob- 
ject G3 10.01 35-hOO. 3892. Th is is an object of pa rticular interest 
as it is the object studied by iKraus et al ] (I2010h via a hish an- 
gular resolution VLTI aperture synthesis image in the ^-band. 
We generate 20 jim images of the set of models that best fit the 
object's SED as found by ^Mo ttram et al.l (1201 ll) . The result is 
shown in Fig. (8] The SED fitting tool preferentially returns mod- 
els with large cavity opening angles (~ 30-50°), but these all 
fail to re produce the observ ed intensity profile. The best fitting 
model of iKraus et al. l (l20Toh. which is also based on SED fitting 
with the model grid of IRobitaille et al.l (l2007l) and the code of 
Whitney et al., provides a better match. However, it still over- 
predicts the flux in the central V from the source (see Fig. (8]). 
This is likely due to the large opening angle of the Kraus et al. 
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Fig. 8. The azimuthally averaged intensity profile of G3 10.0 135 
(upper circles) compared to a PSF standard (lo wer circles) and 
the mo dels that fit the SED of this object from iMottram et alJ 
(I2OIII da shed lines). The so lid line is the radial profile of the 
model of iKraus et al. l (l20Tol) . The error bars represent the rms 
within a given annuli and thus represent an upper limit on the 
uncertainty in the flux distribution. 

model (40°) generating a MIR emission region which is too ex- 
tended. The VISIR observations clearly demonstrate that, for the 
assumed intermediate inclination, a narrow opening angle is re- 
quired. We conclude that the spatially resolved VLT- VISIR im- 
ages at 20 Jim provide valuable constraints on the nature of the 
envelopes surrounding young massive stars. 



7. Conclusions 

In this paper we present difl'raction limited, MIR imaging of 
a sample of massive young stellar objects drawn from the 
RMS survey. By comparing the spatially resolved images of the 
MYSOs to 2D radiative transfer models, we constrain the struc- 
ture of the envelopes that surround them. Here we list the salient 
results. 

- We spatially resolve the MIR (19.5 or 24.5 jim) emission of 
14 MYSOs. The remainder of the MYSOs observed (6) are 
unresolved at the current resolution (-0.6'')- 

- The model of the MYSO W33A that we have developed pre- 
viously can, in most cases, reproduce the images and SEDs 
of the MYSOs. This is not the case for objects associated 
with a Hii region. That the model can reproduce the infrared 
emission of a sample of MYSOs suggests that the circum- 
stellar environments of MYSOs are relatively uniform. The 
large envelope masses of the models are consistent with the 
identification of MYSOs as objects in their accretion phase. 

- It is found that the extent of the MIR emission of the spatially 
resolved MYSOs can generally be attributed to warm dust 
located in the walls of outflow cavities. We suggest that the 
relatively shallow intensity profile gradients found by previ- 
ous MIR diff'raction limited imaging of MYSOs is indicative 
of outflow cavities (rather than alternative explanations such 
as rotational flattening). 

- We show that the varied morphology observed can be at- 
tributed to outflow cavities seen at a variety of inclinations. 
This emphasises that outflows are an ubiquitous feature of 
massive star formation. 
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Appendix A: Notes on individual objects 

Here we discuss the images and modelling results for each re- 
solved MYSO. 



A: G263.7759-00.4281 

The 20 Jim morphology of G263.7759-00.4281 (IRAS 08448- 
4343) is reminiscent of the cavities of a bipolar outflow seen 
close to edge on. The direction of the elongation seen in the 
MIR image is consistent with that seen in the NIR via 2MASS 
images. The radio luminosity of the object is consistent with a 
jet rather than an Hii region (Hoare et al. 2007). This supports 
the notion that the bulk of the MIR emission arises in the walls 
of cavities evacuated by an outflow. Furthermore, the extension 
in the NW/SE direction is aligned with the H2 jet detected by 
iGiannini et al.l (l2005l) , confirming that this object drives an out- 
flow and the MIR morphology is associated with outflow activ- 
ity. 

The SED of this object can be reproduced with the selected 
model with an inclination of / = 87, which is consistent with 
the notion that this object is seen close to edge on. The inten- 
sity profile of the model is not entirely consistent with the data. 
However, the axis- symmetric model used cannot account for the 
slight asymmetry of the image. Indeed, the model can approxi- 
mately recreate the extension of the southern part of the observed 
bipolar structure but not the northern half (see Fig. [71 panel A). 
The model cannot reproduce asymmetry at these wavelengths in 
an edge on configuration. Therefore, we surmise that the model 
reproduces the data as well as can be expected. 

B:G265.1 438+01. 4548 

This object exhibits a fairly symmetric morphology in the MIR. 
If the generic model featuring outflows is correct, this would 
suggest that this object is seen at a relatively low inclination. 
Modelling of this o bject's CO bandhead emission is consisten t 
with this scenario (iBik & Thil [20041: IWheelwright et al.ll2010h . 
The SED and intensity profile are relatively well reproduced by 
the generic model with an inclination of / = 32, which is consis- 
tent with the hypothesis that this object is viewed at a low incli- 
nation. The model slightly under-predicts the flux at distances 
greater than ~2'' from the centre of the intensity distribution. 
However, examination of the image reveals that this may not be 



related to the source as the core of the emission is concentrated 
within T\ Therefore, this slight discrepancy is neglected. 

C: 0268.3957-00.4842 

G268. 3957-00.4842 exhibits a relatively symmetric morphol- 
ogy in the MIR. As before, it is suggested that this is the result 
of a relatively low inclination. Indeed, the object's SED and in- 
tensity profile are reproduced by a model with an inclination of 
/ = 30, consistent with the hypothesis that this object is observed 
at a low inclination. 

D:G269.1 586-01. 1383 

The 20 yum image of G269. 15 86-0 1.1 3 83 incorporates an addi- 
tional, localised source of flux approximately 1'' to the NE. Both 
sources exhibit an extended, slightly cometary morphology rem- 
iniscent of a Hii region. However, radio observatio ns indicate 
that w hile the northern source is a Hii region (Urq uhart et al.l 
[2007h . the southern source is likely to be a YSO. This is con- 
firmed via low resolution NIR spectra. The spectrum of the 
southern source exhibits H2 line emission, suggesting that it 
drives an outflow. 

The complex appearance of this object is difficult to recreate 
with the 2D, axis- symmetric codes of Whitney et al.. In princi- 
ple, a bipolar cometary morphology could be recreated by warm 
dust emission from the walls of cavities carved by a bipolar out- 
flow with a large opening angle. However, the image contains no 
hint of a bipolar structure. Therefore, it might be expected that 
the model cannot reproduce the intensity profile of this source, 
as is found to be the case. We note that this object exhibits a no- 
tably diff'erent morphology to the rest of the sample, which may 
indicate it represents a diff'erent phase of massive star formation 
than the other MYSOs. This is partially substantiated by the fact 
that this object is associated with a Hii region. 

E:G31 0.01 35+00.3892 

G3 10.0 135-^00.3892 (IRAS 13481-6124) was observed with the 
VLTI and AMBER by Kraus et al. (2010). These authors recon- 
structed images of this object in the ^-band and detected a disc 
seen under moderate inclination. In addition jKraus et al ] (l20TQh 
discovered signs of an outflow orientated in the NE/SW direc- 
tion, perpendicularly to the disc. The 20 jim morphology of this 
object is slightly extended in the NE/SW direction, i.e. along the 
direction of the outflow. This is consistent with the notion that 
the majority of the MIR emission traces warm dust in the walls 
of outflow cavities. 

The observed SED and intensity profile are reproduced rel- 
atively well by a model with an inclinati on of / = 3 2 , clos e 
to the intermediate inclination derived by iKraus et al ] (120101) . 
The model under-predicts the flux at distances greater than 
-1.5'' from the peak of the intensity distribution. However, ex- 
amination of the image reveals that at this distance the morphol- 
ogy is not symmetric. The axis- symmetric model cannot repro- 
duce this morphology. Therefore, it is surmised that the model 
reproduces the data as well as possible and this slight discrep- 
ancy is neglected. 

F:G31 8.0489+00.0854 

The 20 fim image of G3 18.0489-^00.0854 exhibits a slightly 
cometary morphology. This object is associated with two Hii re- 
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gions de tected via radio obse rvations, both within several arc- 
seconds (lUrquhart et al.ll200"7l) . The near/far ambiguity over the 
distance to this source has yet to be resolved. As a result, it is im- 
possible to unambiguously model its intensity profile and SED. 
Therefore, we do not attempt to reproduce the observations of 
this object. 

G:G332.2941 +02.2799 

G332.2941-h02.2799 (IRAS 16019-4903) exhibits a signifi- 
cantly extended morphology in the NE/SW direction. This is 
also seen in the NI R in 2MASS images and at 10 jim (see 
iMottram et al .112007b . This object is known to be associated with 
several bipolar outflows and the direction of the most prominent 
outflow is ali gned with the extension seen in the infrared images 
(iHenning et al. 2000). As a result, the observed morphology is 
consistent with the notion that the MIR flux largely arises from 
warm dust in outflow cavity walls. 

Since these observations were taken, it has been deter- 
mined that the lumin osity of G332.2941-h02.2799is 844 L© (see 
iMottram et al .1120 111) . This places the spectral type of this ob- 
ject at approximately B2/B3. Therefore, this object is an inter- 
mediate mass YSO rather than a massive YSO and should be 
excluded from the sample. Nonetheless, we attempt to recreate 
the observations with the generic model developed for the sam- 
ple of MYSOs. The model cannot fully recreate the extension of 
the intensity profile observed, even with a high inclination and 
stellar luminosity. This may imply that the model developed for 
MYSOs is not directly applicable to intermediate mass YSOs. 
Alternatively, the strong outflow may cause shock heating in the 
cavities, which is not included in the model. 

H: G332.9868-00.4871 

G332. 9868-00.4871 exhibits a fairly symmetric morphology. 
As before, it is suggested that this is the result of a relatively 
low inclination. The object's SED is well reproduced by a model 
with an inclination of / = 15°, as is the object's intensity pro- 
file. The model intensity profile deviates from that observed at 
a distance of ~1.75''from the centre. However, this is where the 
object's flux distribution falls to the level of the background and 
thus this discrepancy is not considered significant. 

1:0339.6221-00.1209 

This object exhibits a slightly cometary morphology. As dis- 
cussed in the case of G269. 1586-01. 1383, this is diflScult to 
recreate at 20 jim with an axis- symmetric code. We note that this 
object is the niqst lum inous in the sample with L = 5.2 x lO^L©. 
IMottram et al.l (l201ll) find that there is not a significant popula- 
tion of radio quiet MYSOs with such a luminosity. Sources with 
luminosities greater than L ~ lO^L© are almost inevitably asso- 
ciated with a Hii region. Therefore, it might be expected that this 
object is associ ated with a Hii region . While this object is a radio 
non-detection (lUrquhart et al.l [20071) . it is possible that a weak 
Hii region has escaped detection due to the large distance to this 
object (13 kpc). This could explain the cometary morphology of 
this source. 

A relatively inclined model (/ = 60°) reproduces the object's 
SED and intensity profile relatively well. The mm flux is under 
predicted, but at the large distance to this object, several sources 
could have been included in the mm beam, resulting in an over- 
estimation of the flux. We note that, although the model recreates 



the observed intensity profile, the model image exhibits bi-polar 
structure while the observed morphology is appears more mono- 
polar. As discussed above, this source may harbour a Hii region, 
which could explain this discrepancy. 

J:G343.5024-00.0145 

The image of G343. 5024-00. 0145 exhibits a complex, notably 
extended morphology. Extended nebulosity towards the SE is 
also visible in 2MASS and Glimpse IRAC images. The ex- 
tension of the 20 yum image is broadly consistent with that 
seen at shorter wavelengths, although these data exhibit a more 
complex morphology . This source was detected at 4.8 GHz 
(lUrquhart et al.ll2007h . and is thus classified as a Hii in the RMS 
catalogue. Therefore, it might be expected that the model de- 
veloped to recreate the observations of radio quiet MYSOs can- 
not reproduce the complex morphology exhibited by this object. 
Indeed, the model that recreates the majority of the object's SED 
fails to match the spatial extension of the 20 jim flux. We suggest 
that this is likely due to some of the observed flux originating in 
dust surrounding a Hii region, rather than the central source. 

K:G343.521 3-00.51 71 

The image of G343.52 13-00.5 171 reveals three localised 
sources. The central source clearly dominates the flux. However, 
one of the secondary sources is within 2'' of the primary source. 
It is apparent that secondary flux will contaminate the intensity 
profile of the primary. Indeed, the image reveals that this source 
is less centrally concentrated than its isolated counterparts. A 
limited range of angles was used in constructing the azimuthally 
averaged intensity profile of this object in an attempt to limit 
contamination, but it is not clear how eff'ective this was. The 
model fails to simultaneously reproduce the broad central maxi- 
mum of the image and the observed SED. We suggest that this is 
the result of contamination of the intensity profile by flux from 
the neighbouring secondary source. 

L: G345.0061+01.7944 

G345. 006 1-1-01.7944 exhibits an extended morphology similar 
to that of G332.2941-h02.2799. This suggests that the MIR im- 
age traces outflow structure seen close to edge-on. Indeed, the 
NIR spectrum of this object exhibits H2 emission indicative of 
outflow activity. Unlike the case of G332.2941-h02.2799, a sat- 
isfactory fit to both the SED and the intensity profile was found. 
This is consistent with the hypothesis that the difl&culties in re- 
producing the observations of G332. 294 l-h02. 2799 were the re- 
sult of applying a model developed for high luminosity MYSOs 
to the case of an intermediate mass YSO. 

M:G349.721 5+00.1 203 

G349.7215-h00.1203 is the second most luminous object in the 
sample with L = 3.1 x lO^L©. The image of this object reveals a 
cometary shaped source of flux to the NW of the primary object. 
Radio obse rvations confirm that the additional source of flux is a 
Hii region (lUrquhart et al.ll2007l) . The model reproduces most of 
the features of the SED and intensity profile of this object. The 
model under-predicts the flux at ~100 yum. However, this could 
be due to the measured flux including a contribution from the 
Hii region. The model intensity profile exhibits a slightly diff'er- 
ent slope than that observed, although it is essentially within the 
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uncertainty. We suggest that this could also be due to the pres- 
ence of the Hii region. 
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Fig. 2. Unresolved sources. The images are scaled logarithmically. North is to the top of the page and East is to the Left. The 
contours typically represent 1, 5, 25 and 75 percent of the peak flux. 
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Fig. 3. Resolved sources. The images are scaled logarithmically. North is to the top of the page and East is to the Left. The contours 
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Fig. 3. continued. 
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Fig. 6. The SEDs of the resolved MYSOs compared to the final model SEDs. 
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Fig. 7. The azimuthally averaged intensity profiles of the resolved MYSOs (upper circles) alongside the model radial profiles 
(squares) and the associated PSF standard (lower circles). The lower limit to the y axis is set to the root-mean- square noise in 
the background of the MYSO images. The upper limit of the x is axis is set to the distance at which the MYSO profiles fall to the 
level of the background noise. The radial profile of G263. 7759-00.4281 is also shown averaged over the upper and lower half of its 
bipolar morphology separately (short and long dashed lines respectively). The error bars represent the rms within a given annuli and 
thus represent an upper limit on the uncertainty in the flux distribution. 
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